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Liouville transformations of Schrodinger equations preserve the scattering amplitudes while chang¬ 
ing the effective potential. We discuss the properties of these gauge transformations and introduce 
a special Liouville gauge which allows one to map the problem of quantum reflection of an atom 
on an attractive Casimir-Polder well into that of reflection on a repulsive wall. We deduce a quan¬ 
titative evaluation of quantum reflection probabilities in terms of the universal probability which 
corresponds to the solution of the V 4 = —C 4 /Z* far-end Casimir-Polder potential. 


I. INTRODUCTION 

Quantum reflection of atoms from the van der Waals 
attraction to a surface has been studied theoretically 
since the early days of quantum mechanics [1, 2]. Though 
the classical motion would be increasingly accelerated 
towards the surface, the quantum matter waves are re¬ 
flected back with a probability that approaches unity at 
low energies, due to the rapid variation of the potential 
close to the surface. 

Quantum reflection was first studied experimentally 
for He and H atoms on liquid helium films [3-5] and more 
recently for ultracold atoms or molecules on solid surfaces 
[6-12]. Meanwhile, many theoretical papers have studied 
various fundamental aspects and applications of quantum 
reflection [13-21]. More recently, it has been noticed that 
quantum reflection could be useful for storing and guid¬ 
ing cold antihydrogen atoms [22-25] and that it should 
play a role in any experiment where antihydrogen atoms 
interact with a matter plate (see [26-28] and references 
therein). 

Paradoxical results appear in the study of quantum re¬ 
flection (QR) from the Casimir-Polder (CP) interaction. 
The probability of quantum reflection not only increases 
when the velocity of the incident atom is decreased, but 
also when the magnitude of the interaction is decreased. 
For example, the probability of quantum reflection is 
larger for atoms falling onto silica bulk than onto metallic 
mirrors [26] and is even larger for nanoporous silica [27]. 
This paradox is qualitatively explained by the fact that 
atoms get closer to the surface for a weaker potential, 
so that the CP potential becomes steeper and quantum 
reflection probability larger. 

In the present paper, we propose a quantitative treat¬ 
ment of these paradoxical behaviors based on Liouville 
transformations. Such transformations are gauge trans¬ 
formations of the Schrodinger equation, which map prob¬ 
lems corresponding to different potentials into one an¬ 
other, while leaving scattering amplitudes invariant. In 
the case of QR on a CP potential studied in this paper. 
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a special Liouville gauge can be introduced to transform 
the potential from an attractive CP well into a repul¬ 
sive wall. The paradoxical features of the initial QR 
problem become intuitive predictions of the transformed 
problem. Furthermore, QR probabilities can then be de¬ 
scribed in terms of the universal solution associated with 
the V 4 = —Ci/z'^ far-end Casimir-Polder potential. 

In § II, we recall the usual treatment of the QR prob¬ 
lem, based on deviations from the semiclassical WKB 
approximation [13]. We present in § III the Liouville 
transformations which transform Schrodinger equations 
into equivalent ones corresponding to different potentials. 
We then introduce in § IV a special choice which maps 
the original problem of QR on a CP well into a more 
intuitively understood problem of reflection on a repul¬ 
sive wall. In § V we study the V 4 = —Cilz'^ potential 
which shows non trivial symmetry properties while being 
representative of the CP interaction in the far-end. We 
finally use these results (§ VI) to give a simple evaluation 
of QR probabilities, in terms of the universal function as¬ 
sociated with this problem and of one scattering length 
parameter depending on the full CP potential. 


II. QUANTUM REFLECTION 


We consider a cold atom of mass m incident with a 
velocity u < 0 parallel to the z—axis upon the CP po¬ 
tential V{z) in the half-space z > 0 above the material 
surface located at z = 0. We consider a plane material 
surface, so that the vertical motion is decoupled from the 
horizontal one. 

The vertical motion is then described by a ID 
Schrodinger equation: 


<I'''{z) + F{z)^iz) =0 , 
F(z) . 


( 1 ) 

( 2 ) 


The primes represent the derivative of a function with 
respect to its argument. The CP potential V{z) is at¬ 
tractive, with characteristic inverse power laws at both 
ends of the z—domain, the cliff-side close to the surface 
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z —> 0 and the far-end z —)■ oo: 


V{z) : 

z- 

To ’ 

^ 3 ( 2 ) 

= -C^/z^ , 

(3) 

V{z) 

- Vdiz) , 

Vdiz] 

III 

1 

5 

(4) 


In (1), E = is the energy associated with the 

motion orthogonal to the plane. It may be the result of 
a free fall from a height h above the surface E = mgh, 
where g is the acceleration of gravity. This semi-classical 
treatment of the free fall is possible when E is much larger 
than the energy Ei of the first gravitational quantum 
state above the surface [28, 29]: 


El 



1.407 peV , 


(5) 


where Ai ~ 2.338 is the absolute value of the first zero 
of the Airy function Ai. In the following, we use Ei as 
the unit for E so that the validity of the semi-classical 
treatment of the free fall above the surface is simply 
E/Ei 1. This condition also corresponds to /i/hi 1 
where hi = Ex/mg is the associated unit for the free fall 
height h, that is the height ~ 13.7/xm of the classical 
turning point for the first quantum state [29]. 

If we were treating also the interaction with the CP 
potential in a semiclassical approach, the function F{z) 
would be seen as the square of the de Broglie wave-vector 
kdB associated with the classical momentum hkdB- 

kdBiz) = \/F{z) , ( 6 ) 


As the CP potential is attractive and the incident en¬ 
ergy positive, E is everywhere positive, so that a classical 
particle undergoes an increasing acceleration towards the 
surface. This behavior is mimicked by the semiclassical 
WKB wave-functions which propagate in the rightward 
and leftward directions (77 = -|-I and —1 respectively): 


'^wkb(^) 
OidB{z) = 




\/kdB{z) 


, (/dBiz) 


(7) 

f kdBiz')dz' . ( 8 ) 

J Zn 


Here, adB is the WKB amplitude and (j)dB the WKB 
phase associated with the classical action integral Hc/dB- 
The arbitrariness of the choice of the reference point zq 
is fixed in the following by the far-end convention: 


lim {4>dB{z) - Kz) = 0 , 

Z—¥(X) 

(9) 

- hm kdB{z) - 
2 —>-oo ri 

(10) 


In a quantum treatment of the interaction with the 
surface [26], QR appears as a consequence of non- 
adiabatic transitions between the counter-propagating 
WKB waves (7). It can be obtained by solving the exact 
Schrodinger equation (1) for a wave-function written, in 
full generality, as a linear combination of !fwKB('^) 
z—dependent coefficients Pri'. 


F^z) — /3r;(^)^WKB(^) > (H) 


where we use an implicit sum rule for repeated indices. 
These coefficients obey coupled first-order differential 
equations [13]: 


P'niz) 


P-r,{z) 


2 kdB{z) 


( 12 ) 


This system of equations can be solved numerically and 
matched to the WKB solutions at both ends of the do¬ 
main [18, 30, 31]. Special care has to be taken on the 
cliff-side where the potential diverges. The matching has 
to use the mathematical solutions of ( 1 ) known for the 
V 3 potential [ 22 , 26], at the price of losing physical un¬ 
derstanding of the problem. 

Matter-waves can be reflected back from the cliff-side 
so that the complete problem depends on the details of 
the physics of the surface, including possible sticking, 
non specular reflection or annihilation for antimatter. In 
this letter, we focus our attention on the one-way prob¬ 
lem where the CP potential is crossed only once and, 
therefore, do not discuss this surface physics problem any 
longer. The numerical solution of (1) leads to reflection 
and transmission amplitudes depending on the incident 
energy E or, equivalently, on the parameter n defined 
in (10). A qualitative criterion for occurrence of QR is 
that the coupling term in ( 12 ) takes significant values, 
which may be stated as a large enough variation of the 
de Broglie wavelength on the length scale fixed by the de 
Broglie wavelength. 

A less approximate discussion of this point can be 
based on the remark that WKB wave-functions obey an 
equation differing from the original one ( 1 ): 

'^WKb(^) + ^WKb(^) = ^WKb(^) ' (1^) 

CtdByz) 

The difference between (1) and (13) is often described in 
terms of the so-called badlands function [26]: 

Qiz) = = -aUz)a'Uz) , (14) 

F{z)adB[z) 

which can also be expressed in terms of the Schwarzian 
derivative {^dB,^} of the WKB phase ^dB: 


p. ^ _ {(l>dB,z} _ {<pdB,z} 

2F(z) 2kix,{z) ■ 

The Schwarzian derivative is defined for f{z) as: 

/f m r(^) 

/'(z) 2/'(z)2- 


(15) 


(16) 


An important property of the function Q{z) is that it 
vanishes not only in the far-end where fcdB goes to a con¬ 
stant, but also at the cliff-side as a consequence of the 
power-law variation (3) of the Casimir-Polder potential. 
It follows that the regions where QR takes place are indi¬ 
cated by significant values of the peaked function Q{z). 
This property has been proven by several examples in 
[26, 27] and it is also illustrated on Eig.l. 
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Figure 1. [Colors online] The blue, green and red curves (from 
the highest to the lowest peak) show the functions Q{z) calcu¬ 
lated for energies E respectively equal to 10®, 10^ and 10®Ei 
and V(z) calculated for an hydrogen atom and a silica bulk. 


The three curves show the functions Q(z) for an hydro¬ 
gen atom above a silica bulk with energies E respectively 
equal to 10®, 10"^ and 10®Ei. The curves on Fig.l have a 
higher and higher peak for lower and lower energies, with 
a peak farther and farther from the cliff. This variation 
of the peak value is well correlated with the QR proba¬ 
bility calculated by solving the Schrodinger equation as 
indicated previously (details in [26]). The QR probabil¬ 
ities R, given in Table I, indeed increase when the peak 
values increases, that is also when the energy decreases. 


E ]Ei] 

1 

10 

10® 

10® 

10® 

10® 

R [%] 

98.5 

95.4 

86.1 

63.2 

28.0 

5.6 


Table I. Quantum reflection probabilities R for hydrogen 
atoms falling on a silica bulk plate, with different incident 
energies E, given in units of Ei (see (5)). The three last 
columns correspond to the plots on Fig.l, while the first one 
is given only as an indication of the trend of QR probability 
for the lowest gravitational quantum state. 


In spite of its effectiveness, the method reminded in 
this section suffers several drawbacks. First QR is a 
scattering process with incident matter waves reflected or 
transmitted when crossing the badlands, but this scatter¬ 
ing problem is poorly defined as the potential diverges at 
the cliff. Second the correlation of the peak value of Q 
with the QR probability is observed, but a quantitative 
interpretation of this correlation is missing. In fact, the 
role of the badlands function is discussed by comparing 
the two different problems associated with (I) and (13), 
and not by studying the equation of physical interest (1). 
All these drawbacks are cured in the sequel of this paper, 
thanks to the introduction of Liouville transformations of 


the Schrodinger equation [32]. 


III. LIOUVILLE TRANSFORMATIONS 


The Schrodinger equation (1) is an example of a 
Sturm-Liouville equation under Liouville normal form 
[33]. What is called the WKB approximation by physi¬ 
cists was in fact introduced by Liouville [34] and Green 
[35] for studying properties of such equations, long before 
Wentzel [36], Kramers [37] and Brillouin [38] used it in 
the context of quantum mechanics. The transformations 
used in the present paper were introduced by Liouville 
in 1837 [34] and we follow here the convention of Olver 
[39-41] for naming them after Liouville (see the notes at 
the end of ch. 6 in [39]). 

Liouville transformations have been used to obtain ap¬ 
proximate solutions [42-44]. They have been also used to 
study second-order differential equations applied to solv¬ 
able Schrodinger equations [45-49]. We want to empha¬ 
size at this point that we use here the Liouville transfor¬ 
mations to describe equivalent scattering problems corre¬ 
sponding to different potentials, with no approximation, 
for potentials not belonging to a class of solvable prob¬ 
lems. 

After these historical remarks, we define the Liou¬ 
ville transformations [40] which correspond to coordinate 
changes with correlated rescalings of the wave-function. 
The coordinate change maps the physical 2 —domain into 
a 5—domain with z{z) a smooth monotonous function 
(z'(z) > 0). Equation (1) for 'E{z) keeps the same form 
for the rescaled wave-function <F{z): 


E{z) = ^/¥{^<I'{z) , (17) 

E''{S) + F{i)W{z) =0 , (18) 


with a transformed function F(z): 


m 


F{z) - ^{z,z} 
z'(z)2 


(19) 


The curly braces denote the Schwarzian derivative (16) 
of the coordinate transformation z{z). 

The composition of two Liouville transformations 
z ^ z and z —>■ £ is a Liouville transformation z —>■ z, and 
the group properties of this composition law is ensured 
by Cayley’s identity for Schwarzian derivatives [40]: 


{z,z} = (z'(z))^ {z,z}-k{z,z}. (20) 


When this identity is applied to inverse transformations 
(z = z), the following relation is obtained: 

0=(z'(z)f {z,z} + {S,z}, (21) 

so that the transformation (19) can also be written: 

F(z) = z'(zTF(z) + i{A5} . (22) 
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The Wronskian of two solutions Wi , 'P 2 of the 
Schrodinger equation is a constant, independent of z and 
antisymmetric in the exchange of the two solutions: 

W {'^1,^2) ='^i{z)^2{z) - Ki.zW2{z) . (23) 

The Liouville transformations preserve this Wronskian: 

W(fi,'f'2) = W(#i,'^2) , (24) 

and this property has important physical consequences, 
as shown in the sequel of this section. 

When If' is a solution of (1), its complex conjugate W* 
is also a solution, and the probability density current j{z) 
is proportional to the Wronskian of if"* and S': 

3 {z) = ^W{^\'P) . (25) 

Zim 

As the probability density p = is time-independent 
in the problem studied in this paper, j{z) is a constant 
j independent of z. It follows from (24) that this con¬ 
stant is invariant under Liouville transformations j = j ■ 
We emphasize that the probability density p is neither 
constant nor preserved by the transformation, as one de¬ 
duces from (17) that p{z) = z'{z)p{z). We also note 
that the WKB functions, which are exact solutions of 
(13) and approximate solutions of (1) at both ends of 
the z—domain, obey the following relations (which will 
be used below): 



We recall now that the scattering amplitudes can 
be written in terms of Wronskians of solutions of the 
Schrodinger equation [50]. To this aim, we note that Q{z) 
goes to 0 at the left and right ends of the z—domain, so 
that we can define exact solutions and if"^ of (1) which 
match the asymptotic WKB waves there: 

^l(^) 3 ^„'^wkb(^) 5 ^r(^) 7 ^ '^wkb(^) • (27) 

z—^O z—^00 

These four solutions are schematized in Fig.2. 

As the WKB waves and iF^kb complex con¬ 

jugates of each other, this is also the case for the exact 
solutions matching them at left or right ends: 

(V'^r=v- ( 28 ) 

A generic solution tf'(z) of (1) can then be decomposed 
over the left or right basis: 

^{z) = a^^^l{z)=a^^l{z) , (29) 

For each decomposition, the amplitudes can be collected 
in column matrices related by a transfer matrix: 



These solutions and associated amplitudes can alter¬ 
natively be defined in terms of outgoing and incoming 
waves with the usual identification: 








+ 


out 





Figure 2. Schematic representation of the four solutions which 
match WKB waves on the left and right sides of the interac¬ 
tion region symbolized by the grey ellipse. 


The out and in amplitudes are related by the unitary 
scattering matrix (55^ =71): 

, (32) 

which can be obtained from the transfer matrix (see [51]): 


5 = 

+-I-+1 

r = 



1 

(1 

-r+\ 

' n 


1 ) 


detr= 1 . 


(33) 

(34) 


As the Wronskian of solutions or iF^ is a constant, it 
can be evaluated in particular in the asymptotic regions 
where these exact solutions reduce to WKB waves. They 
therefore obey the same relations as in (26): 

W ((>F(])* , iF(]') = W ((iF^)* , <) = (35) 

The information on the scattering is then contained in 
the Wronskians involving solutions at the left and right 
ends: 

w(('^lT,<) =2 (TT)^'. (36) 

Since the matrices T and S are expressed in terms of 
Wronskians by using (36) and (33), it follows from (24) 
that they are invariant under Liouville transformations 
(T = 'T and S = S). In particular, the reflection and 
transmission amplitudes r = and t = SZ defined 
for waves incoming from the far-end amplitudes are pre¬ 
served (r = f and t = t) and can be calculated equiva¬ 
lently after any Liouville transformations. 

It is worth stressing again that these gauge transfor¬ 
mations relate equivalent scattering problems to one an¬ 
other, while not necessarily making the resolution sim¬ 
pler. In specific cases, for example the model studied in 
§ V, they may lead to non trivial symmetry properties. In 
the general case, we show in the next section that a spe¬ 
cial gauge choice brings satisfactory answers to all points 
raised at the end of § 11. 



5 


IV. SPECIAL GAUGE CHOICE 


In this section, we choose a special Liouville gauge 
which shows interesting properties. Precisely, we choose 
a coordinate z proportional to the WKB phase (/)dB which 
maps the initial problem of QR on an attractive well 
into a different problem of reflection on a repulsive wall. 
This special gauge choice brings satisfactory answers to 
all questions raised above. In particular, it leads to a per¬ 
fectly well-defined scattering problem with no interaction 
in the asymptotic states, and it also allows to understand 
the variation of the QR probability in a rigorous as well 
as intuitive manner. 

The special gauge choice is fixed by the following defi¬ 
nition of the coordinate z and associated quantities iden¬ 
tified by boldfacing: 

Z = ^ ijf(z) = , (37) 

F(z) = 2 =^^l- Q{z)) , (38) 

where we have noticed that {z, z} = {4>dB,z} and then 
used the definition (15) of Q{z). The scale constant k is 
arbitrary at this point, but will be fixed soon. Equation 
(38) can be rewritten in terms of energy and potential: 

F(z) = E - V(z) , E = , V(z) = x^Qiz) . (39) 


As Q{z) goes to zero at both ends of the physical domain 
z G ]0, oo[, the interaction potential V tends to 0 at both 
ends of the transformed domain z g] — oo,oo[. It thus 
corresponds to a well-defined scattering problem with no 
interaction in the asymptotic input and output states. 

Using the expression (14) of Q(z), a positivity property 
can be demonstrated for the integral of V: 

/ oo ^OO 

V(z)dz = -X j adYi{z)ad^{z)dz 

-OO J 0 

pOO 

= xl (adBW)^dz>0, (40) 

where the integrated term, which should appear in the 
integration by parts between the first and second lines, 
vanishes at left and right ends (adBOdB —t 0 for z —0 or 
z —>■ oo). Whereas the initial potential V was everywhere 
negative, the shape of the transformed potential V is 
mostly a repulsive wall, even though V may be negative 
in some parts of the z—domain. We will see now that the 
transformed equation has classical turning points where 
F = 0 or E = V, for not too large values of the original 
energy B. 

Before going further in the discussion of this point, we 
fix the scale constant x to be determined by the wavevec- 
tor K and the length scale £ associated with the far-end 
tail of the CP potential: 


X = VkI , 


y/2rriE 




(41) 


This choice will lead to functions V (z) having nearly 
identical peak shapes for different energies E, at least 
for not too large values of E. In fact, these functions 
reproduce a universal function V 4 (z) when the initial 
potential V{z) matches the form of the far-end tail Vd^z) 
of the CP potential. This model is studied in § V and 
it is shown there that the universal function V 4 (z) has a 
peak value at z = ^ where: 



The plots on Fig.3 show E and V for CP potentials 
V calculated between a hydrogen atom and a silica bulk 
[26] and 3 incident energies E respectively equal to lOEi, 
lO^Ei and lO^Ei. The transformed energies E are rep¬ 
resented as the 3 horizontal lines and the transformed 
potentials V(z) as the 3 curves. With E always posi¬ 
tive and V(z) mostly positive, a logarithmic scale is used 
along the vertical axis, in order to make some details 
more apparent. 



Figure 3. [Colors online] The plots represent the constants E 
(horizontal lines) and the functions V(z) (curves) calculated 
for different scattering problems, corresponding to the same 
CP potential V (z) between an hydrogen atom and a silica 
bulk and energies E = 10, 10^, 10®Ei respectively for the 
blue, green and red curves (from the lowest to the highest 
value of E, or from the lowest to the highest value of V in the 
left-hand part of the plot). The black (dashed) curve is the 
universal function V 4 (z) calculated for a 14 model. 

We see on Fig.3 that the peaks for the functions V are 
nearly the same for the different initial energies. This 
is due to the fact that, for the parameters chosen here, 
these peaks correspond to distances z ^ C such that the 
exact CP potential V{z) is close to its far-end tail V 4 (z). 
The deviations appearing on the plots correspond to val¬ 
ues of z closer to the cliff-side, where V^z) is indeed a 
poor representation of V (z). As V is nearly the same for 
the different problems whereas E = it follows 


K = 
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that classical turning points appear in the transformed 
problem for not too large energies E (x k'^. For the plots 
drawn on Fig.3 turning points appear for E = 10, lO^Ei, 
but not for E = lO^Ei. 

The existence of classical turning points in the trans¬ 
formed problem is in striking contrast with the initial 
problem of QR on an attractive well, which did not show 
turning points. This initial QR problem has been trans¬ 
formed into the more intuitive problem of ordinary reflec¬ 
tion on a repulsive wall, with exactly identical scattering 
amplitudes. The fact that the QR probability goes to 
unity when k —> 0 is now understood as an immediate 
consequence of the increasing reflection expected for a 
particle with a lower and lower energy E coming onto a 
repulsive wall with a more or less fixed peak value. 

In a similar manner, we can understand the depen¬ 
dence of QR probabilities on the absolute magnitude of 
the CP potential. To do so we consider hydrogen atoms 
falling onto a perfect mirror, a silicon bulk or a silica 
bulk, which give rise respectively to weaker and weaker 
CP interaction [26]. Fig.4 shows the constants E and the 
functions V(z) for a fixed energy E = lO^Ei ~ 1.4 neV. 
The potentials correspond to values for the far-end tails, 
C4 and which decrease from perfect mirror to silicon 
to silica. As on Fig.3, the transformed potentials V have 
similar peak shapes, which tend to align on the universal 
curve calculated for a pure V4 potential and shown as 
the dashed curve. In contrast, the transformed energies 
E = nt decrease with £, which immediately explains why 
the QR probability increases [26]. 



Figure 4. [Colors online] The plots represent the constants E 
(horizontal lines) and the functions V(z) (curves) calculated 
for different scattering problems, corresponding to a fixed en¬ 
ergy E = lO^Ei and the CP potentials V(z) for an hydrogen 
atom above a perfect mirror, a silicon bulk and a silica bulk 
(respectively blue, green and red from the highest to the low¬ 
est value of E, or from the lowest to the highest value of V 
in the left-hand part of the plot). The dashed (black) curve 
is the universal function V 4 (z). 


We note that a similar discussion has been given in 
[32] to explain the results obtained in [27] for hydrogen 
atoms above nanoporous silica with different porosities. 


V. SYMMETRY OF THE V 4 MODEL 


In this section, we discuss the model potential 
¥ 4 ( 2 ) = —C^jz'^ which is representative of the CP inter¬ 
action in the far-end. Furthermore, this model is interest¬ 
ing in itself because it obeys a symmetry which enforces 
non trivial properties. 

For the V 4 model, the WKB wave-vector has the simple 
form: 


kdB{z) = a / k 2 + 


£2 

^4 


(43) 


This leads to a non trivial symmetry property for the 
Liouville transformation corresponding to inversion: 


C" 


(44) 


which maps the physical domain z € [0,00] into an in¬ 
verted domain z € [—00,0], while exchanging the roles of 
the cliff-side and far-end. 

The inversion (44) is an homographic function, so that 
its Schwarzian derivative {z, z} vanishes. If z is another 
map, chosen arbitrarily, Cayley’s identity (20) leads to 
{z, z} = (z'(z))^ {z, z}. When z = one deduces that 
the badlands function defined as in (15) for the original 
and inverted coordinates are identical: 


Q{z) = 


{(l>dB,z} ^ WdB,z} 

2fcj3(z) 


= Q{i 


(45) 


with ^dB{z) = (j)dB{z), fcdB = <?^dB(2^)> ^dB = ^dB(5)- 
The badlands function may be written explicitly: 


Q{z) = 




2«2 




2/12 


(k2z2 -I- I 2 ) (fi;2z2 -I- I 2 ) 


= Q{z) , (46) 


and it reaches its peak value at z = ^, that is also 
z = ( = —(. This peak value scales as the inverse of 
= k£. When multiplied by the latter value (see (41)), 
it leads to the universal function: 



5 

8cosh^(2M) 


u = \n^ , (47) 


with a peak value | and a relation between z and u still 
to be discussed. 

The WKB phase and the coordinate z also obey sym¬ 
metry properties under the inversion: 


z = 


= [ \/2 cosh(2'u') du' , 

Juo 

pU _ 

z* -I- / 1/2 cosh(2u') du' , (48) 
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with z* the value corresponding to the inversion center: 

z* = z(C) = ^r(|)\ (49) 

Eqs (47) and (48) constitute an explicit parametric rep¬ 
resentation of the universal function V 4 (z) proving that 
it is symmetrical with respect to the inversion u —^ —u, 
that is also z — z* —>■ — (z — z*). Another representa¬ 
tion in terms of hypergeometric functions is given in the 
Appendix B, where other homogeneous forms of the po¬ 
tential V{z) are also considered. 

The QR probability calculated for the 14 model [52, 
53], denoted i ?4 in the following and plotted on Fig. 5, is 
a universal function of the dimensionless parameter It 
can be calculated by numerically solving the Schrbdinger 
equation for the potential Vi{z) or V 4 (z). Alternatively, 
it can be obtained by the analytical method summed up 
in the Appendix A. 



Figure 5. Quantum reflection probability R 4 , calculated for 
the 14 model and shown as a function of the dimensionless 
parameter 


VI. DISCUSSION OF QR PROBABILITIES 

We now present the values obtained for QR proba¬ 
bilities, and compare the exact results for the full CP 
potential with those obtained for the V4 model. 

We first recall that the QR probability goes from unity 
at K —>■ 0 to zero at k —>■ 00 . Its departure from unity at 
low energies is described by a scattering length a defined 
by the general relation: 

r{K) ~ ~ (1 ~ 2jKa) . (50) 

K —>-0 

The scattering length is a complex number, the imag¬ 
inary part of which determines the quantum reflection 
probability: 

R{k) = \r{K)\'^ ~ 1 — Anb , b = —?sa. (51) 


Table II gives i and b for an hydrogen atom above a 
perfect mirror, a silicon bulk and a silica bulk, as ob¬ 
tained from the full calculations in [26]. The table shows 
that the equality b = £ typical of the 14 model [23] is no 
longer true for the full CP potential, with large variations 
in particular for the case of silica bulks. 


mirror 

perfect 

silicon 

silica 

b [ 00 ] 

543.0 

435.2 

272.6 

£ [ao] 

520.1 

429.8 

321.3 


Table II. Comparison of the values of b and £ for an hydrogen 
atom above a perfect mirror, a silicon bulk and a silica bulk, 
given in atomic units ag — 52.9 pm. 


We show on Fig.6 the calculated QR probabilities i? 
as a function of the dimensionless parameter nb for the 
scattering problems corresponding to Table II. The full 
curves represents the values calculated for perfect mir¬ 
rors, silicon and silica bulks [26]. They are compared 
to the dashed curve which corresponds to the universal 
function R 4 (with b = £ in this case) calculated for the 
pure 14 model. Using the value calculated for b for dif¬ 
ferent bulks of matter, it turns out that the exact QR 
probabilities R are close to the expression R 4 evaluated 
for the same value of nb. The agreement is excellent at 
low energies because the peaks of V on which QR occurs 
correspond in these cases to values of z in the far-end 
tail of the CP potential. Differences between the curves 
R and Ra appear at large enough energies, because dis¬ 
tances z closer to the cliff-side play a significant role. 



Figure 6 . [Colors online] Log-log plot of the quantum reflec¬ 
tion probability R shown as a function of the dimensionless 
parameter nb, corresponding respectively to a perfectly re¬ 
flecting mirrors (blue), silicon (green) and silica (red) bulks 
(from the lowest to the highest curve at the right-hand side of 
the frame) and compared to R 4 calculated for the V 4 model 
(black dotted curve). 

The same information is given in Table III with 












precise numerical values of the QR probabilities R 
calculated for the same scattering problems at en¬ 
ergy E = lO^Ei ~ 1.407 neV. This corresponds to 
K = 8.237X 10® m“^ that is also 4.359 x 10“^ Oq ^ (atomic 
units with oq — 52.9 pm). The comparison with 
obtained for the calculation of the V 4 model shows good 
agreement in accordance with the fact that QR occurs 
in these cases in the far-end tail of the CP potential. 
It is worth stressing that the agreement would be much 
poorer when comparing R to R 4 {k£). 


mirror 

conductor 

silicon 

silica 

nb 

0.237 

0.190 

0.119 

R [%] 

41.8 

49.2 

63.2 

Rdinb) [%] 

41.9 

49.0 

63.1 


Table III. Quantum reflection probabilities R for hydrogen 
atoms falling on a perfectly conducting, a silicon and a silica 
bulk plate, with incident energy E = lO^Ei, compared to the 
corresponding value of the universal function i? 4 (K 6 ). The 
values of the b and t are given in Table II. 


In this paper, the problem of QR of an atom on a 
Casimir-Polder attractive well has been mapped into an 
equivalent problem of reflection on a wall through a Li- 
ouville transformation. This gauge transformation of the 
Schrodinger equation relates exactly equivalent quantum 
scattering processes which correspond to different semi- 
classical pictures. It produces a new interpretation of 
the main features of quantum reflection and explains in a 
clear manner the paradoxical features of the initial prob¬ 
lem. It also allows quantitative evaluation of QR proba¬ 
bilities which can be obtained from the universal function 
corresponding to the pure V 4 model. 

Acknowledgements - Thanks are due for insightful 
discussions to M.-T. Jaekel, V.V. Nesvizhevsky, A. Yu. 
Voronin, and the GEAR and GRANIT collaborations. 


Appendix A: QR probability for the I4 model 

In this appendix, we recall the analytical method which 
can be used to to solve Schrodinger’s equation for the V 4 
model [52, 53]. The derivation presented here follows the 
work of [54] and uses results in [40, 55]. 

We first perform a Liouville transformation: 

z^z(z)=ln^, E(z)^E(S) = ^. (AI) 

C V ^ 

With the new variables, the Schrodinger equation for the 
V 4 model becomes a modified Mathieu equation: 

'^"{z) + (—a + 2q cosh(2z)) W{z) = 0 , (A2) 

where a = | while q = x = y/nl is the only remaining 
parameter. A pair of solutions to this equation can be 


written as series involving products of Bessel functions: 

00 

n——oo 

(A3) 

Here r is a complex parameter yet to be determined, 
known as the Mathieu characteristic exponent. The co¬ 
efficients An '’ obey the following recurrence relation: 

((t -I- 2n)^ - a) A^l'^ + q (^A^J+i + ^n-ij = 0 • (A-4) 

The infinite determinant associated with this sys¬ 
tem of equations must be zero for a non trivial so¬ 
lution to exist. This singles out a value of r, 
which can be obtained by following the procedure 
detailed in [54] or using the Mathematica function 
MathieuCharacteristicExponent[a,q] . With the re¬ 
currence relation we can write the ratios An'’ / and 
as continued fractions. These ratios go to 
0 when jn| increases so that we can truncate the contin¬ 
ued fractions to obtain numerical values for An'’ (with 

A’ = 1 ). 

As a result of the invariance of equation (A 2 ) under 
parity z —>■ —z (which is the symmetry discussed in § V), 
Ei±){-z) are also solutions, and one can show that: 

#(=*=) ( 5 ) = (j = lnf—(A5) 

^ ^ ^ ^ <f'(+)( 0 ) 

Using known results for the Bessel functions: 


Jn{x) ~ \ —cos (, 

x^oo V TTX \ 


X -T-7 


^ ; 'Aj(O) — dnfi 5 


we deduce the asymptotic behaviors: 

¥^\z) 



As we are looking for the reflection and transmission 
amplitudes r and t for a wave coming from the far-end, 
we search the solution {z) of Schrodinger equation 
(A2) which has the asymptotic behaviors: 


m-{z) cs 


tz 

0 v^' 


exp ( —i ( — 


m^{z) ~ 


(A7) 


We have used the asymptotic forms of (j)dB{z) 


(f'dBiz) ~ 2xZt, -, (l)dB{z) ~ KZ 

2—>-0 z 2—>-00 


(AS) 
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Matching the asymptotic forms (A6-A7), we obtain the 
reflection and transmission amplitudes: 


r = 


—i 


sinh(CT) 

sinh(cr + iitt) ’ 


sin(7rT)e^*^^* 
sinh(tT + iitt) 


(A9) 


It has been checked that this analytical method gives the 
same results as a direct integration of the Schrddinger 
equation [26]. The resulting QR probability R 4 = jrp is 
drawn on Fig. 5. 


We define the special Liouville gauge z = (j)dB{z)/>c as 
in (37) and obtain E = and V = k^Q^z) as in (39). 
Up to know, we did not fix the scale factor as in (41), as 
£ = £4 does plays not play any role for the potential U„ 
studied in this Appendix. 

As the maximum value of Q scales as l/(/cC„)^, we 
choose the scale factor as: 

Xn = nCn , (B7) 


Appendix B: Other homogeneous potentials 


In this appendix, we consider the case of homogeneous 
potentials Vn{z) = —Cn/z" with n > 2. This includes 
the already discussed case n = 4, as well as the case n = 3 
which corresponds to the Van der Waals zone close to the 
surface and the case n = 5 which corresponds to the far- 
end of a slab mirror [26]. 

We start by introducing two relevant length scales: 


= -UW .(Bl) 

They generalize the definitions of C and i for n = 4. 
Qn and in measure respectively the distance at which 
E =\Vn\ and the strength of the potential. 

The WKB wavevector and phase are thus read: 

kAB{z) = n\ I + \ , x=^, (B2) 

V CC" Cn 

(t>dB{z) = nCn [ \/1 + ^ dx' , (B3) 

Jxo * ^ 

where Xq is chosen to enforce (9). For n > 2, we deduce: 


<(>dB = 


nXK(^n 
n — 2 

\ \ z n n 


(B4) 



where F is the hypergeometric function dehned as in [40]. 
The particular case n = 4 gives an alternative expression 
for the expressions in § V. 

The badlands function: 


nx” ^ 4 — n + 4(1 + n)x” 
^ {<uY 16(l + x")^ 


(B5) 


is a peaked function reaching its maximum at: 


I h'n? — 3n — 8 + — 6 n^ — 13n^) 


4(n2 + 3n + 2) 


(B6) 


which generalizes the definition (41). This leads to uni¬ 
versal functions V„ which do not depend on any other 
parameter than n. The functions V„(z) are drawn on 
Fig.7 for the cases n = 3,4, 5, as functions of the coordi¬ 
nate z = (j)^B{z)/K using the convention (B7). 



z 


Figure 7. Plot of the universal functions V„(z) for n = 3, 4, 5 
(smallest to tallest). 


With the same convention, the integrals (40) of V„ 
over the real axis are real numbers depending only on 
n and they can be expressed in terms of the Gamma 
function: 



In particular, the case n = 4 corresponds to: 


. _ 5F(5/4)2 


0.772531 . 


(B8) 


(B9) 
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